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Conversion of natural habitat into agricultural landscape has been identified as one of the major drivers 
of habitat loss. Human-modified ecosystems, such as agricultural land, have gained significant attention 
in terms of the conservation of their native biodiversity. We studied the effectiveness of organic agroe- 
cosystems in conserving odonate diversity by comparing organic terrace rice cultivation with a nearby 
natural forest system with streams in Sikkim, Eastern Himalaya, India. We sampled adult odonates using 
a transect count method (laying six permanent transects) covering two villages (Lingmoo in South Sikkim 
and Dzongu in North Sikkim), making a total of 48 transect counts. A total of 881 individual odonates 
representing 31 species under two suborders (16 Anisoptera and 15 Zygoptera) and seven families were 
recorded during this study. Of these, 20 species representing three families and 18 species represent- 
ing seven families were observed in terrace rice cultivation and the natural forest system, respectively. 
Beta diversity estimates showed that the community composition of the odonates differed qualitatively 
(incidence measure) and quantitatively (abundance measure) between the two land use types. Turnover 
component (abundance balance in case of abundance based beta diversity measure) had higher contribu- 
tion in the overall beta diversity, suggesting that one assemblage of species is being replaced by another 
due to environmental sorting. The variation in community composition between the two habitats was 
statistically significant. Our results suggest that organic wetland habitats are important for conservation 
of odonates and associated biodiversity (especially herpetofauna) in the Himalaya and require urgent 
conservation attention. 
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Introduction 


Habitat loss and degradation have resulted in significant decline in biological diversity over the 
last decade (Brooks et al., 2002; Primack, 2014). Conversion of natural habitats into agricultural 
landscapes has been identified as one of the major drivers of habitat loss (Newbold et al., 2015) 
and a high percentage of red-listed species is known to be threatened by agricultural intensifica- 
tion (IUCN, 2015; Norris, 2008). Agricultural intensification is predicted to increase in the future 
to meet the food demand for the ever-increasing human population, resulting in higher numbers 
of threatened species and loss of natural ecosystems (Godfray et al., 2010). Environment friendly 
farming practices, which often involve traditionally and organically managed agroecosystems, 
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are some of the promising approaches that can maintain food security while safeguarding bio- 
diversity. Such farming practices support higher species richness and abundance of organisms 
than conventional farming (Czech, 2002; Katayama et al., 2019; Lichtenberg et al., 2017; Tuck 
et al., 2014). The effectiveness of environment friendly farming practices often depends on use 
of organic approaches and availability of surrounding landscape as a source habitat for range of 
biodiversity (Baba, Kusumoto, & Tanaka, 2019; Tscharntke, Klein, & Kruess, 2005). 

Habitat alteration due to agriculture intensification has been identified as one of the key drivers 
of local extinction of arthropods (Chung, Eggleton, Speight, Hammond, & Chey, 2000; Claus- 
nitzer et al., 2009; Lawton et al., 1998). Intensive expansion of rice fields is considered as the 
main reason for loss of habitat of many wetland species (Russi et al., 2013). Among aquatic 
arthropods, odonates represent a group which is severely affected by intensification of rice culti- 
vation (Kadoya, Suda, & Washitani, 2009). Odonates depend on a large number of environmental 
factors and their species assemblages are determined by different characteristics of the habitat, 
including water quality (Corbet, 1999; Gibbons, Reed, & Chew, 2002; Scher & Thiery, 2005). 
Alteration of natural habitat causes the replacement of forest specialists by generalist groups of 
odonates (Clausnitzer, 2003; Samways, 2003; Seidu, Danquah, AyineNsor, Amaning Kwarteng, 
& Lancaster, 2017). Studies from southern Brazil (Dalzochio, Baldin, Stenert, & Maltchik, 2016), 
north-western Italy (Giuliano & Bogliani, 2019) and Japan (Baba et al., 2019) have indicated the 
positive effect of environmental friendly practices of rice cultivation on odonate communities. 

In India, rice cultivation is one of the primary agricultural practices and almost 55% of its over- 
all population consumes rice (Jadhao & Khurad, 2011). Rice is grown in 43.86 million hectares 
in India, making it an important semi-natural landscape (Ministry of Agriculture and Farmers 
Welfare, 2019). This landscape acts as a habitat for odonates to survive and reproduce, where 
natural habitats are unavailable. A total of 12 species of odonates were recorded from the rice 
fields of Madurai, Tamil Nadu (Kandibane, Raguraman, & Ganapathy, 2005), 21 species from 
Tiruvallur district of Tamil Nadu (Anbalagan, Paulraj, & Ignacimuthu, 2013), 13 species from 
Howrah district of West Bengal (Dwari and Mondal, 2017) and 14 species from Anand in Gujarat 
(Rohmare, Rathod, & Parasharya, 2016). The diversity of odonates, however, was recorded to be 
lower in agricultural landscapes than in other natural systems (Kulkarni and Subramanian, 2013). 
Conventional methods of intensive agriculture are practised in almost every state in India. For 
example, 7596 of farmers in Tamil Nadu use agrochemicals for their rice cultivation (Chitra, 
Muraleedharan, Swaminathan, & Veeraraghavan, 2006). The use of agrochemicals, namely fer- 
tilizers and pesticides, had a significant negative effect on the population of odonates in rice 
fields in Karaikal, Puducherry, India (Soniyagandhi & Kumar, 2014), indicating that intensive 
rice cultivation practices pose a serious threat to odonate diversity in agricultural landscapes. 

The agroecosystem in Sikkim, a small state of India, comprises of an agroforestry system 
with forest habitat surrounding the agricultural landscape (Sharma, Chettri, Gurung, & Shakya, 
2007). Sikkim was recently declared as an organically farming state in India (Bhutia, 2015). The 
agroforestry system of Sikkim was found to be effective in conserving high butterfly diversity, 
including species of conservation concern (Chettri, Sharma, Dewan, & Acharya, 2018, Sharma 
et al., 2020). The rice agroecosystem in Sikkim comprises of about 10.67 thousand hectares of 
cultivated land (Sharma & Acharya, 2013). These human-modified systems in the region can be 
crucial and serve as an alternative habitat for biodiversity, especially water dependent species 
including odonates. 

The Odonata fauna of Himalaya is highly diverse, with 257 species in 112 genera and 18 
families (Subramanian and Babu, 2017). The high diversity of Odonata is reported from Eastern 
Himalaya, especially in Sikkim-Darjeeling and Arunachal Pradesh (Mitra, 2003, 2004; Allen, 
Molur, & Daniel, 2010). From the Sikkim-Darjeeling region alone, a total of 185 species of 
odonates including 13 endemics are reported (Allen et al., 2010; Payra & Bhutia, 2017; Subra- 
manian & Babu, 2017). Additionally, the region is recognized as one of the global freshwater 
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biodiversity hotspots (Allen, Molur, & Daniel, 2010). Traditionally managed farming systems 
have been proved to be an alternative refuge for various taxa, both invertebrate and verte- 
brate (Kurniawan, Yanuwiadi, Priambodo, Maulidi, & Kurnianto, 2016; Sharma et al., 2020). 
The terraced farming system in Sikkim also acts as refuge for various species of herpetofauna. 
Snake species such as Indo-Chinese rat snake (Ptyas korros), Himalayan keelback (Herpetoreas 
platyceps), checkered keelback (Xenochrophis piscator) and some frog species of the genus 
Fejervarya and Polypedates live in sympatry with odonates in paddy fields. These animals are 
considered as farmers’ friends, as they control various pests such as rodents and insects. A recent 
study from lowland Nepal highlighted that frogs are potential biological control agents and thus 
an integral part of the rice paddy system (Khatiwada et al., 2016). 

Earlier studies on odonates in the region focused mainly on taxonomic diversity and distri- 
bution. However, no systematic studies on association of odonates with agricultural systems 
have been carried out in Sikkim, where agricultural practices are fully organic. Therefore, this 
study was conducted to document the diversity of odonates in the rice agroecosystem of Sikkim 
Himalaya and to understand the effectiveness of such agroecosystems in conserving odonate 
diversity. Specifically, we compare the alpha and beta diversity of odonates between organic 
terrace rice cultivation (TRC) and nearby natural streams in the forests. 


Methodology 


Study area 


Sikkim is a small hilly state of India with an area of 7096 km? and is a part of the larger Himalaya 
biodiversity hotspot. The steep rise in elevation in this region causes adverse changes in climatic 
profile, resulting in huge variation in floral and faunal assemblages (Acharya, Sanders, Vijayan, 
& Chettri, 2011). Two villages in Sikkim, Lingmoo and Dzongu, were selected for the present 
study (Figure 1). Lingmoo is a village in the South District of Sikkim, located at an elevation of 
1100 m. The vegetation consists of sub-tropical mixed broad-leaved hill forest. Dzongu village 
is located in the North District of Sikkim and has an elevation ranging from 800 m to 6100 m. 
We considered elevations between 800 and 1000 m for this study because TRC studied here was 
available only within this range. The region supports the different vegetation types; however, 
the typical vegetation of our study sites was sub-tropical broad-leaved forest. The agricultural 
ecosystem of these regions typically consists of cultivated land interspersed between patches of 
natural forest. The details of these land use systems are given below. 


Terrace rice cultivation (TRC) This system mostly consists of cultivation of rice in artificially 
created terraced lands (Figure 2). Unlike in mainland plains of India, rice in Sikkim is grown 
in water-fed terraced steep terrain. Indigenous varieties of rice, which only grow in the ambient 
climatic conditions of the region, make TRC a unique agricultural area and important ecosystem. 
Due to its extensive irrigation and water drainage system, TRC forms an important artificial 
habitat for odonates and some species of amphibians and reptiles. 


Natural forest system with streams (NFS) The landscape of the study area consists of patches of 
privately owned forests, community managed forest or continuation of protected forest around 
the vicinity of agricultural land (although protected forests were not considered in the present 
study). The local community chiefly depends on these forests for fuelwood, fodder, herbal 
medicine and also for some wild edibles. Forest around the elevation of 900-1200 m (the range 
covered during this study) mostly comprises of Alnus nepalensis, Schima wallichi, Castanopsis 
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Figure 1. Map showing the study sites of Odonates in Lingmoo (South Sikkim) and Dzongu (North Sikkim), Eastern 
Himalaya. 





(C) (D) 


Figure 2. Photographs showing two habitat types: terrace rice cultivation (A, B) and natural forest with streams (C, D) 
considered for the study in Sikkim, Eastern Himalaya. 
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Table 1. Details of transects laid for the study of Odonates in Sikkim, Eastern Himalaya, India. 


Transect Village Elevation (m) Latitude (°N) Longitude (°E) Habitat type 





Tl Lingmoo 965 27.33590 88.46672 NFS 
T2 Lingmoo 1133 27.34113 88.46270 TRC 
T3 Lingmoo 1043 27.34037 88.46495 TRC 
T4 Dzongu 816 27.52120 88.53145 NFS 
T5 Dzongu 956 27.50283 88.51978 TRC 
T6 Dzongu 998 27.52695 88.51335 NFS 


Abbreviations: NFS, natural forest system with streams; TRC, terrace rice cultivation. 


indica, Engelhardtia spicata, Macaranga indica and Terminalia myriocarpa, as dominant veg- 
etation. Lantana camara, Rubus niveus and Boehmeria sp. mainly dominates the undergrowth 
vegetation. Water bodies such as streams and rivulets present in this ecosystem are mostly undis- 
turbed and free from anthropogenic disturbances (Figure 2). The natural water bodies provide 
habitat to aquatic and amphibious organisms that depend on a moist forest environment. 


Odonata sampling 


We sampled adult odonates using a transect count method following the line of least resistance. 
Three permanent line transects of 200 m length and 5m width were laid in both the TRC and 
nearby streams in the NFS, making six transects in total (Table 1, Figure 1). Transects were laid 
along artificial rivulets, water channel and paddy pools in the TRC system. Natural forest areas 
with streams were selected as reference sites for comparison of odonates between disturbed and 
reference habitats. Transects in NFS were laid along the stream line with forest cover. A slow and 
steady pace was maintained during the transect walk and all the odonates found within 5m band 
on both sides of the transect were recorded through visual observation. Sampling was carried out 
between 09:00 and 14:00 hours during calm and sunny weather days because most odonates are 
active and the probability of encounter was high during this time. While most individuals were 
identified visually without capturing, species which could not be identified easily (e.g. Pseuda- 
grion sp. and Calicnemia sp.) were captured with the help of a sweep net, and released after 
taking notes on essential morphological characters. The odonates which remained unidentified 
in the field were photographed and later identified with the help of the official portal of odonates 
of India (Joshi, Dawn, Roy, & Kunte, 2019) and available literature (Fraser, 1933, 1934, 1936; 
Mitra, 2003; Subramanian, 2005; Subramanian & Babu, 2017). Sampling was conducted during 
the months of July 2018 to April 2019 covering all major seasons. A total of 25 and 23 transect 
count were made in the TRC and NFS respectively. No collection of the samples for preservation 
was made during the study. 


Data analysis 


Community parameters such as species richness, abundance, and alpha and beta diversity indices 
were calculated for odonates of the study area. Species richness was considered as the total num- 
ber of species encountered during the survey and abundance as the total number of individuals 
of the odonates. Non parametric estimators were also computed to estimate the richness using 
EstimateS software (version 9.1.0) (Colwell, 2013). Among the set of estimators Chaol was 
selected owing to its high precision in estimating richness (Hortal, Borges, & Gaspar, 2006). We 
then generated a species accumulation curve using this estimator and observed richness in order 
to ensure the completeness of the sampling. The Shannon—Wiener diversity index (H^) was cal- 
culated as the component of alpha diversity (Pielou, 1969; Shannon & Weaver, 1949) using the 
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formulae given below: 
Shannon-Wiener diversity index (H^) = — Xpjlnp; (1) 
where p; — proportion of the total sample belonging to the ith species; In — natural logarithm. 
Evenness (J ) = H'/Hgax (2) 


where H4, = In$; $ = number of species; H' = Shannon—Weiner diversity index. 

To compare the alpha diversity and its associated parameters between TRC and NFS, we fol- 
lowed two distinct approaches. First, we applied the Mann-Whitney U test (the data assumed 
non-parametric distribution) to test the significance in variation between the community param- 
eters of TRC and NFS. Second, in order to estimate the magnitude of difference in diversity 
of the two communities, the effective number of species (ENS), also known as Hill's number, 
was estimated. Effective number of species is estimated as the exponent of the Shannon-Weiner 
diversity index and represents the number of equally abundant species required to produce a 
particular value of diversity. 

Beta diversity between TRC and NFS was calculated based on incidence and abundance 
matrix and partitioned into its respective components. To explain beta diversity in terms of 
incidence, Sorensen dissimilarity index (sor) was calculated for total beta diversity, Simpson 
dissimilarity (Psm) as turnover component and nestedness resultant (Byes) as its nestedness com- 
ponent, following Baselga (2010). For abundance-based dissimilarity measurement, total beta 
diversity was calculated as Bray-Curtis dissimilarity (dgc) and then partitioned into balanced 
variation in abundances (dpc.o4), which complements species turnover and abundance gradi- 
ents (dgc-gra) as a complement of species nestedness (Baselga, 2013). Permutational multivariate 
analysis of variance (PERMANOVA) and non-metric multidimensional scaling (NMDS) was 
performed based on Bray—Curtis dissimilarity to understand the variation between the TRC and 
NFS Odonata community assemblages. All statistical analysis was performed in R software ver- 
sion 3.4.3 using vegan (Oksanen et al., 2016) and betapart (Baselga et al., 2018) packages. We 
then summarized the differences in community by comparing the community assemblages at 
family level between the TRC and NFS, because the differences are more pronounced at family 
level as suggested by Sigutova, Sipos, and Dolny (2019). 


Results 


Diversity and abundance 


A total of 881 individual odonates representing 31 species in two suborders and seven families 
were recorded during this study (Appendix 1; Figures 3, 4). Suborders Anisoptera and Zygoptera 
represented 16 and 15 species, respectively. TRC consisted of 20 species representing three 
families and NFS had 18 species representing seven families (Appendix 1; Table 2). Estimated 
species richness (Chaol) was higher in the TRC (20 + 0.31) compared to NFS (19.49 + 2.59). 
The values of observed species richness and Chaol were almost similar for TRC, indicating 
completeness of sampling. Chaol estimated relatively higher values of species richness than 
the observed richness of the odonates in NFS, indicating that a few additional species may 
be recorded on further sampling. Species accumulation curves for both TRC and NFS almost 
reached an asymptote, indicating completeness of the sampling effort (Figure 5). Other com- 
munity parameters such as Shannon- Wiener diversity index, evenness and effective number of 
species (ENS) differed between the two communities (Table 2). However, the difference in the 
community parameters between the study systems was not statistically significant. The ENS 
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Ischnura aurora Rhinocypha unimaculata 


Figure 3. Photographs of some of the damselflies (Zygoptera) recorded during the study in Sikkim, Eastern Himalaya 





Orthetrum sabina 





Palpopleura sexmaculata Orthetrum glaucum Trithemis aurora 


Figure 4. Photographs of some of the dragonflies (Anisoptera) recorded during the study in Sikkim, Eastern Himalaya 
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Table 2. Community parameters of odonates of terrace rice cultivation (TRC) and natural forest system with streams 
(NFS) in Sikkim, Eastern Himalaya. 











SI no Community parameter TRC NFS Total 
1 Species richness 20 18 31 
2 Total abundance 639 242 881 
3 Shannon- Wiener diversity 2.63 2.47 2.55 
4 Evenness 0.88 0.84 0.86 
5 Chaol + SD 20 + 0.31 19.49 + 2.59 - 
6 Effective number of species 13.96 11.91 
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Figure 5. Species accumulation curve in (A) terrace rice cultivation (TRC) and (B) natural forest system with streams 
(NFS) in Sikkim, Eastern Himalaya. S(est)-observed species; Chaol-non parametric species richness estimator. 


value of the TRC and NFS were 13.96 and 11.91, respectively, indicating that the TRC was 
slightly more diverse than the NFS. A higher number of individuals (n = 639) was recorded 
from the TRC as compared to the NFS (n = 242). In the TRC, Palpopleura sexmaculata was 
the most abundant species, constituting about 15.02% of the total odonates observed, followed 
by Orthetrum glaucum (12.37%) and Orthetrum triangulare (11.90%). Similarly, Anisopleura 
lestoides was the most abundant species, comprising 24.80% of the total odonates recorded in 
the NFS, followed by Calicnemia eximia (9.91%) and Orthetrum glaucum (9.09%). 
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Table 3. Beta diversity between terrace rice cultivation and natural forest system 
with streams based on incidence based and abundance based dissimilarity matrices 
in Sikkim, Eastern Himalaya. 







Incidence based Abundance based 
Indices dissimilarity dissimilarity 
Turnover/abundance balance 0.588 (Bs) 0.61 (dgc-bal) 
Nestedness/abundance gradient 0.05 (nes) 0.17 (dac-gar) 
Total beta diversity 0.63 (B son) 0.78 (dgc) 
Unidentified. species. 1 
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Figure 6. Non-metric multidimensional scaling plot showing clustering of terrace rice cultivation (TRC)/natural forest 
system with streams (NFS) habitat transects and species assemblage (based on Bray- Curtis dissimilarity) 


Beta diversity pattern 


Both incidence based (B sog = 0.63) and abundance based matrices (dgc = 0.78) of beta diver- 
sity measurement showed high diversity between the two habitats (Table 3). Partitioning of 
beta diversity showed that the turnover component contributed most to beta diversity. Non- 
metric multidimensional scaling plot (stress — 0.092) showed distinct clusters of two habitat 
types depicting distinct species assemblages (Figure 6). Among the 31 observed species, seven 
species were common to two habitat types, whereas 13 species were unique to TRC and 11 were 
exclusive to NFS (Appendix 1). Species of the genera Calicnemia and Coeliccia made distinct 
clusters of forest species. Similarly, generalist species such as Orthetrum triangulare, Orthetrum 
glaucum, Crocothemis servilia, Diplacodes trivialis, Ischnura aurora and Palpopleura sexmac- 
ulata formed a distinct cluster inhabiting the TRC system. PERMANOVA analysis showed 
statistically significant variation in odonate community between the two habitats (p < 0.001) 
(Table 4). 
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Table 4. PERMANOVA results for terrace rice cultivation and natural forest system 
with streams in Sikkim, Eastern Himalaya based on Bray- Curtis distances. 


d.f. SS MS Pseudo F Rab 
Habitat 1 1.7358 1.73578 3. 0.24334% ** 
Residuals 23 5.3974 0.23467 0.75666 
Total 24 7.1332 1 


Abbreviations: SS, sum of squares; MS, mean sum of squares. 
* Significance value based on 999 permutation. 
^ Significance level: ** p < 0.001. 


Family wise pattern 


Odonates representing seven families, Libellulidae, Gomphidae, Euphaeidae, Platycnemididae, 
Coenagrionidae, Chlorocyphidae and Platystictidae, were recorded during the study. Odonates 
belonging to all seven families were observed in NFS. Platycnemididae comprised of 38.8996 
of total species recorded in this habitat, represented by seven species, followed by Libelluli- 
dae (22.2296) with four species, and Chlorocyphidae and Coenagrionidae (11.1196) with two 
species each. Euphaeidae, Gomphidae, and Platystictidae comprised 5.5696 of species, with one 
species each. In terms of relative abundance, Platycnemididae was the most abundant (37.2096 
of total individuals), followed by Libellulidae (28.9296), Euphaeidae (24.8096), Coenargrionidae 
(7.03%), Platystictidae (0.82%), Chlorocyphidae (0.82%) and Gomphidae (0.41%). 

Of the seven families, only three (Libellulidae, Coenagrionidae and Euphaeidae) were 
recorded in the TRC. Libellulidae was the most dominant family in the TRC (75%), repre- 
sented by 15 species, followed by Coenagrionidae (20%) with four species and Euphaeidae (5%), 
represented by only one species. A similar pattern was observed in relative abundance of the fam- 
ilies, with Libellulidae being the most abundant family (79.03%), followed by Coenagrionidae 
(20.03%) and Euphaeidae (0.94%). 


Discussion 


This is the first systematic study on odonate communities of the agricultural landscape of 
Sikkim, an organically farming state of India. The total species recorded during the present study 
(31 species) represent about 40.2% of the total number of species recorded from the Sikkim 
Himalayan region (Mitra, 2003, 2004; Payra & Bhutia, 2017). Non-parametric estimators esti- 
mated slightly higher number of species than observed species richness, indicating the chance of 
encountering a few more species with additional sampling. However, the asymptotic pattern of 
species accumulation curve indicates completeness of the sampling. 

Alpha diversity and abundance were slightly higher in the TRC compared to NFS. Other stud- 
ies have also reported higher abundance and richness of odonates in agricultural land compared 
to forest areas (Hofhansl & Schneeweihs, 2008; Seidu et al., 2017). However, the differences in 
alpha diversity parameters between the two communities were statistically non-significant in our 
study. This suggests that the two communities are characterized by almost equitable patterns of 
alpha diversity, with TRC slightly more diverse than the NFS. The diversity of odonates within 
the TRC was maintained mainly due to the influx of species with high adaptability to open 
habitat conditions and their tolerance to moderate levels of habitat disturbances, even though 
forest species were being lost. Hence, the inference based only on alpha diversity measures 
(species richness count and its derivative) may not be adequate to compare the two communities, 
because the real differences were in the distinct species compositions of the two communities. We 
observed high beta diversity between the habitats which are mostly dominated by substitution 
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components reflecting distinct community assemblages. Similar results have been obtained by 
a study in Mexico which shows that the land use conversion affects the species composition 
of adult odonates but not the species richness (Rocha-Ortega, Rodrigurz, & Cordoba-Aguilar, 
2018). 

The beta diversity estimates showed that the community composition of the odonates dif- 
fered qualitatively (incidence measure) and quantitatively (abundance measure) between the 
two land use types. The turnover component (abundance balance in case of abundance based 
beta diversity measure) had a higher contribution in the overall beta diversity, suggesting that 
one species assemblage is being replaced by another due to environmental sorting (Baselga, 
2010). This high percentage of dissimilarity indicates the uniqueness of the two habitat types 
as shown in the NMDS plot (Figure 6). NFS is a natural habitat with streams, and harbours a 
specialist community with less tolerance to environment changes, whereas the TRC habitat is 
widespread, supporting rich micro fauna and common species with frequent occurrence of indi- 
viduals that can survive in various environmental conditions (Elanchezhyan, Sowmiya, Agilesh, 
& Venkatesh, 2017). TRC in this region is characterized by landscape modification to support rice 
cultivation. Altered habitats tend to alter species assemblages of odonates and favour generalist 
species which can tolerate wide range of habitat type (Clausnitzer, 2003; Stewart & Samways, 
1998). Species richness was slightly higher in TRC compared to NFS but high turnover high- 
lights the importance of both ecosystems in maintenance and conservation of odonates and 
associated biodiversity. Similar alpha diversity and high beta diversity between the systems high- 
lights that human-modified ecosystems such as TRC are equally effective as forest ecosystem in 
conservation and maintenance of odonate diversity. 

In the TRC system, the odonate community was mostly dominated by generalist species. 
Libellulidae of the suborder Anisoptera was found to be most dominant family and species such 
as Palpopleura sexmaculata, Orthetrum glaucum and Orthetrum triangulare of the same family 
were most abundant. These species are reported to have wide distributional range and can breed 
in wide array of microhabitats, but often prefer marshy lentic habitats (Dow, 2009a, 2009b, 
2010a). Orthetrum triangulare in particular is known to be tolerant to disturbances and can be 
considered as an indicator of habitat alteration (Dow, 2010a). Additionally, most of the species 
under family Libellulidae are swift fliers, which requires open areas for active thermoregulation 
to compensate for the energy lost during their flight (Corbet & May, 2008; Seidu et al., 2017). 
The dominance of Libellulidae in the rice cultivation ecosystem has been reported by many other 
studies as well (Dwari & Mondal, 2017; Kulkarni & Subramanian, 2013; Rohmare et al., 2016). 
Open habitats such as TRC offer more favourable thermal microhabitats than forest, where the 
canopy of trees offers a cooler environment. Similarly, the Coenagrionidae family of the sub- 
order Zygoptera were represented by species such /Jschnura aurora and Aciagrion occidentale, 
which are associated with open and disturbed habitats (Dow, Rowe, & Marinov, 2013; Mitra, 
2010; Sigutova et al., 2019). 

NFS harboured the highest diversity of Zygoptera among the recorded odonate families. Most 
zygopterans requires vegetative substrate for oviposition and hence are favoured in habitats with 
good vegetation cover (Seidu et al., 2017). In NFS, Platycnemididae was the most dominant 
family both in terms of species richness and abundance. The species of this family, such as Cal- 
icnemia sp. and Coeliccia renifera, were restricted to forest habitat because of their dependency 
on the forest for breeding (Dow, 2010b, 2010c). Similarly, the two families of Zygoptera, namely 
Platystictidae and Chlorocyphidae, were rare and strictly confined to forest habitat. Anisopleura 
lestoides, belonging to family Euphaeidae, was found to be particularly dominant in the NFS but 
was found to occur in the TRC (mostly in the margins with some vegetation cover) in less abun- 
dance. Since Euphaeidae can tolerate habitat modification to some extent, it may completely 
disappear only with extensive alteration of natural habitat (Subramanian, Ali, & Ramchandra, 
2008). 
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In summary, our study assessed odonate diversity in the TRC and NFS in the Sikkim Himalaya. 
We found significant differences in community assemblages of odonates in these two land use 
types. Our results indicate that nature-friendly farming practices with remnant forest patches 
in the rice cultivation landscape are able to sustain the diversity of forest specialist groups of 
odonates. In contrast to intensive agriculture practices around the world, the traditional and 
organic farming practice in the Sikkim Himalayan region ensures the preservation of remnant 
forest patches in its surroundings. Considering the intensification of agriculture in tropical coun- 
tries, the role of remnant patches in conserving the biodiversity is well appreciated across the 
globe (Hansen, Scheele, Driscoll, & Lindenmayer, 2019). We are, however, uncertain on how 
different environmental, spatial or biotic variables affect odonates in this landscape, since we did 
not evaluate these factors during the study. It has been reported that gradients of forest degrada- 
tion, area of forest patch and level of disturbances causes shift in odonate communities (Dolny, 
Barta, Lhota, & Drozd, 2012; Luke et al., 2017; Sigutova et al., 2019). Subsequently, it has been 
reported that insecticide usage, water quality, ditch size and structure of bank vegetation in the 
rice agroecosystem significantly affect the odonate community (Baba et al., 2019; Giuliano & 
Bogliani, 2019). Hence, it is important in the future to understand these processes and to identify 
key components that shape community assemblage of the odonates in this agricultural landscape. 


Conclusion 


The study highlighted the effect of habitat changes on odonate assemblages in organic agri- 
cultural landscapes of Sikkim Himalaya. The TRC system had higher species richness and 
abundance which was mostly dominated by the generalist species. On the other hand, NFS had 
the richer family level diversity with most of the species specialized to forest habitat. The rice 
cultivated areas thus represent an example of an altered habitat. However, the adjoining forest 
patches also form an integral part of the agricultural landscape and compensate for the loss of 
habitat for odonates, and serve as a refuge for the forest specialist species. 

It is to be noted that the lowland areas of Sikkim have only 6 km? area under protected area 
network (PA) and further expansion of PAs is not possible in these human-dominated landscapes. 
The remnant forest patches retained within traditional cultivated systems in Sikkim Himalaya can 
be an important alternative to the PAs in conservation of biodiversity. Further studies document- 
ing various taxa such as herpetofauna in this agricultural ecosystem, and the ecosystem services 
associated with them, may be crucial for planning conservation strategies. Studies on structural 
and species composition, effective area and disturbance gradients of the forest patches and their 
interactions would be crucial for future conservation planning. While our results suggest that 
organic agricultural farms are important for conservation of odonates and associated biodiver- 
sity, extensive studies involving more traditional farms across the state (or in the other parts of 
Eastern Himalaya) would provide us with a better understanding of the biodiversity conservation 
potential of these unique ecosystems in the Himalaya. 
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Appendix 1: List of odonates observed during the study in Sikkim, Eastern Himalaya 


Abbreviations: TRC, terrace rice cultivation; NFS, natural forest system with streams (NFS); LC, least concern; DD, data 
deficient; VU, vulnerable (LC, DD and VU as categories as per IUCN 2019). 4- sign represents present and — represent 
































absent. 

Habitat 
SI no Scientific name TRC NFS Total abundance IUCN status 
Suborder Anisoptera 
Family Libellulidae 
1 Orthetrum luzonicum (Brauer, 1868) + - 30 LC 
2 Orthetrum triangulare (Selys, 1878) 97 LC 
3 Diplacodes trivialis (Rambur, 1842) 70 LC 
4 Orthetrum glaucum (Brauer, 1865) 101 LC 
5 Neurothemis intermedia (Rambur, 1842) + - 12 LC 
6 Pantala flavescens (Fabricius, 1798) + - 35 LC 
J Orthetrum sabina (Drury, 1770) + - 8 LC 
8 Crocothemis servilia (Drury, 1770) + - 39 LC 
9 Rhodothemis rufa (Rambur, 1842) + - 9 LC 
10 Lyriothemis sp. + - 3 LC 
11 Palpopleura sexmaculata (Fabricius,1787) + - 96 LC 
12 Orthetrum pruinosum (Burmeister, 1839) + + 51 LC 
13 Trithemis aurora (Burmeister, 1839) + - 7 LC 
14 Orthetrum taeniolatum (Schneider, 1845) + - 14 LC 
15 Unidentified species 1 + - 3 
Family Gomphidae 
16 Anisogomphus bivittatus (Selys, 1854) - + 1 DD 
Sub order Zygoptera 
Family Platycnemididae 
17 Coeliccia didyma (Selys, 1863) - + 10 LC 
18 Coeliccia renifera (Selys, 1886) - + 15 LC 
19 Calicnemia nipalica (Kimmins, 1958) - -+ 10 VU 
20 Calicnemia miniata (Selys, 1886) - + 17 LC 
21 Calicnemia eximia (Selys, 1863) - + 24 LC 
22 Calicnemia pulverulans (Selys, 1886) - + 3 LC 
23 Calicnemia erythromelas (Selys, 1891) - + 11 LC 
Family Coenagrionidae 
24 Pseudagrion sp. 48 LC 
25 Ischnura aurora (Brauer, 1865) + - 44 LC 
26 Ceriagrion fallax (Ris, 1914) + - 7 LC 
27 Aciagrion occidentale (Laidlaw, 1919) 46 LC 
Family Chlorocyphidae 
28 Aristocypha cuneata (Selys, 1853) - + 1 DD 
29 Rhinocypha unimaculata (Selys, 1853) - + 1 LC 
Family Platystictidae 
30 Drepanosticta carmichaeli (Laidlaw, 1915) - + 2 LC 
Family Euphaeidae 
31 Anisopleura lestoides (Selys, 1853) + + 66 LC 


